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The results of measurements and numerical simulation of charge carrier distribution and energy 
states in strained quantum wells In^jGai-^As/GaAs (0.06 ^ x ^ 0.29) by C-V-profiling are pre- 
sented. Precise values of conduction band offsets for these pseudomorphic QWs have been obtained 
by means of self-consistent solution of Schrodinger and Poisson equations and following fitting to 
experimental data. For the conduction band ofi'sets in strained lna;Gai_j, As/GaAs - QWs the ex- 
pression /S.Ec{x) — 0.814a; — 0.21a;'^ has been obtained. 

PACS numbers: 73.21. Fg, 81.07.St 
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INTRODUCTION 



Since the development of semiconductor heterostruc- 
tures the determination of energy band discontinuities of 
various semiconductor pairs has been a very important 
task. Energy band offsets dominantly control the elec- 
tronic states in heterostructures and, hence, the output 
parameters of semiconductor devices. The importance 
of getting true values of band offsets as well as the dif- 
ficulties in obtaining and, even more, in interpreting the 
relevant data have been attracting attention for the last 
30 years. R. Dingle was one of the first who reported in 
1974-75 ^'^ the value of band offsets for isoperiodic het- 
erosystem (Al-Ga)As/GaAs ("Dingle" rule 85:15). Then 
H. Kroemerj^i^ G. Duggan^ and Yu et aL^ comprehen- 
sively reviewed the understanding of band offsets before 
1991 and provided an overview of the methods commonly 
used in experimental band offset determination, mostly 
optical at that time. At the same time, the authors;^ and 
others showed that a low sensitivity of the optical transi- 
tion energies to the band offsets made its determinations 
rather confusing. Up to now a great number of papers 
has been published on this subject (see bibliography in 
recent comprehensive review^). So far, however, as was 
pointed out in the review, among the ternary alloys used 
in quantum electronics, only the AlGaAs/GaAs system 
has generally accepted values of band offsets. 

For one of the most important used heteropairs - 
Ina;Gai_a:As/GaAs - as yet no clear picture about the 
dependence of band offsets on alloy composition has been 
obtained, despite the very intensive investigations in last 
years. The data collected by P. Bhattacharya^^ show a 
great scatter of the values of relative conduction band 
offset ASc between 35% and 85% for x < 0.35. Above 
mentioned reviewS reports relative conduction band off- 
sets for the Ina;Gai_a;As/GaAs - system in the range 57- 
90% and recommends as a rule of thumb AEc[eV] = x 
for X < 0.5. They conclude that no detailed study 



has yet been carried out on InGaAs-based heterojunc- 
tions. Recent publications on this subjecbii*iSiiii4 only 
present partial results for different compositions, more 
or less agreed with "recommended" in Ref. 9. Theoretic 
calculations^^ give the valence band offset for the end 
combination InAs/GaAs AEy = 0.06 eV, which is in 
serious disagreement with experimental data. 

One important device application of the heterosystem 
InGaAs/GaAs is high power laser diodes with strained 
quantum wells^ In these structures thin quantum-size 
layers of InGaAs grow pseudomorphically, i.e. having the 
lattice constant of the underlying GaAs-layer in the plane 
of the heterojunction. The elastic energy, accumulated 
due to crystal cell distortion, causes the band structure of 
the thin InGaAs layer to be modified, altering partic- 
ularly its energy gap. Hence, in strained InGaAs/GaAs 
quantum wells one should expect another band offsets 
than in heterostructures with thick layers of the solid so- 
lution. In cases between pseudomorphic growth and full 
strain relaxation (occurs in thick layers) the band offset 
in InGaAs/GaAs will have, obviously, some intermedi- 
ate values. This fact explains, we suppose, the variety of 
data found in literature. 

Numerical fitting of C-V-curvesi2ii2iS£ by means of 
self-consistent solution of Schrodinger and Poisson equa- 
tions is one of the most promising approaches to measure 
band offsets of quantum well structuresiii*2iiS^ This ap- 
proach correctly takes the quantization of carriers in a 
quantum well into consideration and yields very accu- 
rate results. However, well-defined heterostructures are 
necessary for this. 

Complex multilayer structures like multiQWs etc. and 
unknown dopant profile or the presence of deep levels add 
sources of uncertainties. Therefore, in order to be sure to 
get precise values for band offsets at hetero junctions sim- 
ple structures with a minimum of unknown parameters 
or parameters to be fitted should be used. 

This work presents accurate data for band off- 
sets in heterostructures with strained pseudomorphic 
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Ina;Gai_xAs/GaAs (0 < a; < 0.3) quantum wells. To 
obtain these values we have carried out a systematic cy- 
cle of C-V-measurements on specially fabricated struc- 
tures. Details of sample preparation and measurements 
are described in Section ^ In Section IIIII the model 
for simulating measured concentration profiles and de- 
riving the values of conduction band offsets based on 
self-consistent numerical solution of the Schrodinger and 
Poisson equations is described. The carrier concentration 
in the quantum well region is calculated on the base of 
a quantum-mechanical approach. Mathematical aspects 
of the computations are presented in Sections IIII Al and 
IIII Bl To increase the accuracy of the numerical calcula- 
tions a non-uniform mesh with the mesh step inside the 
quantum well 10 times smaller than in the other regions 
has been used. Finally, in Section llVI we present the re- 
sults of numerical fitting of experimentally measured C- 
V curves. The dependence of conduction band offset for 
strained pseudomorphically grown In2;Gai_2,As/GaAs - 
QWs has been obtained as a function of quantum well 
composition in the range 0.06 ^ x ^ 0.29. 



II. SAMPLE STRUCTURE AND 
MEASUREMENT PROCEDURE 

A special set of high quality samples with a simplified 
structure (Fig. containing Ina;Gai^xAs/GaAs quan- 
tum wells of different width {w = 6.0 — 9.5 nm) and com- 
position {x = 0.065 — 0.29) was grown on n+-GaAs sub- 
strates by MOVPE at deposition temperatures of 650°C 
and 770°C. The GaAs cladding layers were uniformly 
doped with Si, except for the QWs themselves and thin (5 
nm) spacer layers on both sides of the quantum well. To 
get the best experimental results and to eliminate possi- 
ble uncertainties in the subsequent numerical fitting, the 
cap GaAs layer was designed to be 300 nm thick and 
have a constant doping level of (6 — 7) * 10^^ cm~^. The 
width and composition of the QWs and cladding layers 
have been determined by high resolution X-ray diffrac- 
tion (HRXRD). All QWs were fully strained without any 
relaxation seen in X-ray area mapsii^ Ag-Schottky bar- 
riers were fabricated on top of the structures and Ohmic 
contacts were formed on the substrate. 

The parameters of the grown structures are listed in 
Table n 

The measurements of capacitance-voltage characteris- 
tics and profiling of majority carriers in the quantum 
wells have been carried out with the help of a computer- 
controlled G-V-profilometer at a testing frequency of 1 
MHz and with an amplitude of the probing signal of 15 
or 50 mV. 

At zero bias the width of the space charge region un- 
der the Schottky-barrier in the samples was less than the 
thickness of the cap GaAs-layer. With increasing reverse 
bias the space charge region was broadened and its bor- 
der crossed the quantum well. The C-V-characteristics 
of all samples clearly exhibit a plateau in the range of 
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FIG. 1: Layer sequence of the grown samples with 
In^jGai-a; As/GaAs - quantum wells. 
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FIG. 2: C-V-characteristics of Ino.225Gao.775As/GaAs QW 
at different temperatures (Sample #307). 



Uicv = 2 — 4.5 V related to discharging carriers in the 
QW. A typical example of 1/C'^ = f{Urev) characteristic 
for sample #307 at different temperatures is shown in 
Fig. 13 

The apparent carrier distribution is derived from a 
measured C-V-curve using the well known formula for 
depletion approximation 



n{d) = 2 
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where e is the dielectric constant (assumed to be equal for 
both the well and the barriers), e is the electron charge, 
A is the area of the Schottky diode. The depletion width 
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of conduction band discontinuities in 




^As/GaAs strained quantum wells grown by MOVPE. 
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Figure 13 shows some examples of the n — d curves 
covering the whole range of QW compositions. The pro- 
files exhibited a clear dependence of amplitude, width, 
and the depth of depletion on the composition and the 
width of the QWs. It is worth to note that beyond the 
regions of accumulation and depletion related to the QW 
the carrier concentration was excellently constant, and 
we used it in the fitting of the experimental profiles to 
the simulated ones on the base of self-consistent solution 
of the Schrodinger and Poisson equations by varying the 
band offset. 



III. MODEL FOR SIMULATING C-V 
PROFILES 

For simulating the C-V-characteristics the Poisson 
equation 
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has been solved, where ^p{z) is the electrostatic potential, 
n{z) is the free carrier concentration, and is the con- 
centration of ionized donors. The boundary conditions 
for Q at the Schottky barrier and in the electroneutral- 
ity region, far away from the QW, are: 



V3(0) = U + ifbi, 



f{oo) = 0. 



(4) 



(5) 



Here U is the applied voltage, and ipbi is the built-in 
potential. 

In addition, the matching conditions for the potential 
at both heterointerfaces have to be fulfilled 
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The indexes "barr" and "well" correspond to the regions 
of GaAs barrier and InGaAs quantum well, respectively. 

The free carrier concentration n{z) in far from the 
QW can be calculated as in the case of a homogeneous 
structure through the Fermi integral 



n{z) = Nc^F, 
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kT 



(7) 



where Nc is the effective density of states in the con- 
duction band, Ep is the Fermi level, T is the temper- 
ature, and k is the Boltzmann constant. In contrast, 
in the vicinity of a quantum well the carrier concentra- 
tion should be calculated by solving Schrodinger's equa- 
tion. The needed spatial distribution of the electro- 
static potential was derived using a procedure of self- 
consistency of the Schrodinger and Poisson equations. 
The essence of the procedure is the sequential (step-by- 
step) solution of Schrodinger and Poisson equations until 
convergenceiSii22i24 As a criterion of convergence we took 
an increment of the potential of less than 10^® V in the 
next iteration. 

Quantum size effects are important only inside the 
quantum well and in its immediate vicinity. So for nu- 
merical solution of the Schrodinger equation we used 
a " quantum box" , a narrow region containing the QW 
(Fig. 21). The optimal width of the quantum box was 
chosen experimentally during simulations, in order to 
achieve a high precision in determination of the quan- 
tized carrier concentration, and, in balance, to reduce the 
computing time. The quantum box width was chosen as 
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FIG. 3: C-V-profiling of In^Gai-a; As/GaAs quantum wells. 
Samples 7^296-#309. a) Common view of apparent carrier 
concentrations for several samples, T=300K; b) Apparent 
concentration peak values as a function of x. 
Note that near x = 0.14 there are 4 samples with different 
width of QW. 



nine times the width W of the corresponding quantum 
well, which was placed in the center of the quantum box. 

According to the boundary conditions, there must be 
nodes of the wave functions at the edges of the quantum 
box. For this reason the data calculated close to the 
quantum box boundaries are dropped. The length of 
this region is no more than 1 W ^ as can be seen from 
Fig. 21 On the other hand, at distances about 2 — 3 
W from the QW the quantization effect is very weak, 
and there we can use the Fermi integral JT)) for deriving 
the free carrier concentration. The coincidence of the 
concentration profiles at this part derived from quantum- 
mechanical approach and from was used as the proof 
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FIG. 4: Scheme of computations. The "quantum box" is 
shaded. 1 — Schrodinger's concentration; 2 — concentration 
derived from 0; 3 — resulting calculated "apparent" con- 
centration (matched the experiment). 



for a true solution. 



A. Solving the Poisson equation 

The Poisson equation has been solved numerically by 
Newton's method relative to the correction term. 

The great difference (tenfold) in the values of electron 
and hole effective masses in GaAs and nearby ternary 
InGaAs alloys makes the Fermi level shift toward the 
bottom of the conduction band. Therefore, a significant 
part of donors (up to 30%) remains non-ionized at room 
temperature, even despite the very low ionization energy 
of Si donors in GaAs (about 5 me\i2^). Because of this, 
the incomplete donor ionization has to be taken into ac- 
count. 

To reduce the computation time, it is desirable to use 
some approximation of the Fermi integral iQ . The sim- 
plest exponential approximation is not applicable here 
because of the close position of the Fermi level to the 
bottom of the conduction band. There is another well 
known approximation for Q by the expression 



n(z) 



exp 



Ec-EF-eip(z) 
kT 



(8) 



that better matches the Fermi integral. The constant C„ 
here usually falls between 0.17 and 0.35>^ 

In order to minimize the approximation error and to 
fulfill the electroneutrality condition on the right-hand 
side of the simulated region (i.e. in the GaAs substrate) 
we used the following procedure: by solving the elec- 
troneutrality equation and using expression Q the Fermi 
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level position is determined at 1^9 = 0. Then equating (UJ 
to © one can derive the current adaptivity constant C„. 
At another Lp the maximum relative error of such approx- 
imation does not exceed 3 x 10~^. 

The electrostatic potential was written as an initial 
approximation ^po{z) and a correction term lS.Lp[z): 



(9) 



To linearize the Poisson equation Q the expression for 
n was decomposed into a Taylor series including linear 
term relative to the correction l^Lp{z). 

Then a finite-difference analog of the Poisson equation 
has been rewritten as a system of linear equations with 
a characteristic three-diagonal form. To get a high pre- 
cision solution in a reasonable time different mesh steps 
were used inside and outside the quantum box. The num- 
ber of points in the mesh was 8000, including about 1500 
in the quantum box. The Gauss method was applied to 
solve the system with some modifications based on obvi- 
ous symmetry of the equations. 

After getting the correction /S.ip{z) a new potential was 
obtained according to @ • 



B. Solving the Schrodinger equation 

The effective mass, one-dimensional Schrodinger equa- 
tion can be written as^ 



d 1 dijji{z) 
2 dz m*{z) dz 



+ V{z)^^{z)^E,ij,{z), (10) 



where Ei are the eigenvalues, ipi are the corresponding 
eigenvectors, m* is the coordinate-dependent electron ef- 
fective mass. V{z) is the effective potential energy: 
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(11) 



AEc is the conduction band offset. 

We used boundary conditions of Neuman's type at the 
ends of the quantum box. 

The finite-difference analog for (|10|) was obtained using 
the three-point formula 



2m* 



+ Vj^P,^, = E,^P,^,, (12) 



where j identifies the point on the one-dimensional mesh, 
and hj is the distance between the mesh nodes (a step of 
the mesh). 

In addition to boundary conditions, at the heterojunc- 
tions the following matching conditions should be main- 
tained between the derivative of wave functions inside 
and outside the quantum well: 
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The Schrodinger equation H10|l was solved numerically 
by the well known "shooting" method with some im- 
provements aimed to reduce the computation time. 

The number of points in the mesh should be enough 
to eliminate the error due to substitution of the deriva- 
tive with the finite-difference approximation (|12|l . We 
compared the results of numerical solution of (|10|l by the 
shooting method with the well known analytical solution 
for a rectangular quantum welliSi It was found that the 
mesh size of about 1500 points yields quite good accu- 
racy with a relative error in eigenvalue determination less 
than 10^"^ for almost all levels. 

After the set of eigenvalues Ei and corresponding 
eigenvectors 4^i{z) had been obtained, the carrier con- 
centration in the region of QW was calculated via local 
density of states from the expressionSLSS 



m*{z)kT^ 



1 + exp 



Ep — Ei 
kT 



(14) 



using the condition of normalizing the wave functions 



|V'.(z)pdz = 1. 



(15) 



Summation in H14|l runs over all subbands. 

The prefactor in the expression H14|) before the square 
of wave function is considered as the number of electrons 
per unit area in the ith - subband. 

In the quantum well region the concentrations of 
both bound and free electrons were calculated from the 
Schrodinger equation. We consider this a more correct 
approach than the simple summation of bound (inside 
QW) and free carriers (above QW) used, for example, 
in Ref. 22, But, due to the finite width of the quantum 
box the continuum of free electron states in this scheme 
of computations is represented by a set of discrete levels 
with energies determined by the size of the quantum box. 

In order to sum up all charge carriers in the quantum 
well region we took into account the 16 lowest energy 
levels. Fig. |31 shows the occupation of the 16 energy lev- 
els at different bias. The carrier concentration in the 
first subband ni was at least ten times greater than in 
the second one, and the concentration in 16th subband 
(and all higher subbands) is below 10"'' of ni and can be 
neglected. 

The results of the computation of energy states in de- 
pendence of applied reverse bias for the sample #300 
[AEc — 175 meV) are shown in Fig. |HIa), and a lineup 
of the conduction band bottom for this structure is de- 
picted in Fig.EIb). As can be seen, a single bound level 
is observed in the structure with an energy of 49 meV 
in equilibrium. (The bottom of the conduction band in 
the electroneutrality region was taken as zero). Starting 
approximately at —2.5 V the space charge regions of the 
Schottky barrier and the QW merge (Fig. [71), and the 
penetrating electric field bends the conduction band bot- 
tom near the QW, forcing the bound level to lift up. At 
U = —5 V the level becomes unbound. 
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FIG. 5: Electron concentration in first 16 energy subb- 
bands (relative units) at different reverse biases. Sample #300 
(ASc = 175 meV), T=300K. 



To calculate the C-V characteristics so called "qua- 
sistatic approach"^* was applied. The capacitance of a 
structure is the first derivative of the total charge. The 
latter can be derived via the flow of electric field across 
the surface according to the Gauss theorem. The spatial 
distribution of electrostatic potential is calculated dur- 
ing solution of the Poisson equation, so one can derive 
the value of the electric field at the surface at different 
applied biases: 



-^surf — 



^bulk 



(16) 



and, hence, build up the capacitance- voltage characteris- 
tic (or restore the apparent concentration profile) using 
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FIG. 6: Results of numerical calculations for sample #300 
{x = 0.23; well width = 7.2 nm; offset AEc = 175 meV): 
a) first 8 energy levels as a function of Urcv', b) tlie bottom of 
conduction band near the QW at different Uiav 



IV. RESULTS OF SIMULATION AND 
DISCUSSION 

As has been established earlier^SSiSS there exists a cer- 
tain discrepancy between the true and "apparent" con- 
centration profiles of free charge carriers near a het- 
erojunction, a quantum well or a quantum dot^^'^^^ 
An apparent profile, obtained in experiment, is more 
smeared in comparison to the true one and has a shift 
in the peak position (see Fig. 0J. The general reason 
for this discrepancy is the indirect and non-equilibrium 
procedure of concentration profile restoration from C-V- 
measurements. Generally, this technique involves differ- 
entiation of the C-V-curve |^ in the approximation of 
fully depleted space charge region and does not take into 
account the problem of Debye smearing. In the case of 



QW profiling, where one expects sub-Debye resolution, 
this standard technique leads to an essential distortion 
of the apparent profile. So, for the goal of adequate 
fitting, during simulations we must accomplish just the 
same procedure of restoring the apparent profile as in real 
experiment and, particularly, the bias voltage increment 
in the simulation must be equal to the voltage step used 
in the experiment. 

The results of fitting for two samples are presented in 
Figs. IHl and 1^1 As can be seen, excellent matching is ob- 
tained. This proves the correctness of the used model. 
One should underline again that due to the high quality 
of the specially fabricated for C-V-measurements samples 
no additional adjustable parameters like an impurity con- 
centration gradient or a charge at the heterojunction had 
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FIG. 7: Simulated concentration profiles of electrons in the 
region of QW at U = —1 Y (solid) and —2.5 V (dashed line). 
Sample #300. 



to be used in the fitting procedure. The only fitting pa- 
rameter was the conduction band offset, AEc- The value 
of majority carrier concentration was taken on the shoul- 
ders of the measured concentration profiles. Parameters 
for Ina;Gai_a;As, needed for the calculations, were taken 
fromiSS*^ Fig. |5| also demonstrates the resolution of our 
fitting. For a medium In-content (x = 0.14) the error 
was less than 10 meV. It was found that the resolution is 
approximately directly proportional to the alloy compo- 
sition of the quantum well. In general, we estimate the 
relative error in the determination of band offsets as less 
than 10% within the measured range of x. 

An interesting example of fitting for the sample with 
the smallest In-content in QW {x = 0.065, sample #298) 
is presented in Fig. ^| Here the apparent peak of en- 
richment in the QW is even smaller than the value of the 
impurity concentration, despite the spatial confinement 
inside the QW. The simulated profile in this case is very 
sensitive to the band offset (the error is about 5 meV), 
however, the fitting is not as good as for other composi- 
tions. For such a weak concentration peak the presence of 
residual impurities in the quantum well and in adjacent 
spacers begins to play an essential role. One also should 
bear in mind for explanation the increased relative value 
of the experimental noise. 

In Table m we summarized the conduction band offsets 
in strained pseudomorphically grown Inj;Gai_a;As/GaAs 
(0.06 ^ a; ^ 0.29) quantum wells obtained in our study. 
Only one bound level was observed in all samples. Its 
depth in equilibrium ([/ = V) is depicted in Fig. [TTl as 
a function of composition. One can see that for composi- 
tions X < 0.25 the level appears above the corresponding 
Fermi level. Nevertheless, the occupation in the subband 
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FIG. 8: Experimental (dots) and fitted (solid) apparent 
concentration profiles of Ino.23Gao.77As/GaAs quantum well. 
(Sample #300, T=300K). Best fit AEc = 175 meV. 
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FIG. 9; Resolution of the fitting. The results for the best 
fit [AEc = 120 meV, solid line) and for AEc = 130 meV 
(dashed line) are presented. Dotted curve — the experimental 
apparent profile for sample #303 {x = 0.145). 
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FIG. 10: Apparent concentration profile of 

In^^Gai-i, As/GaAs quantum well with low In-content 
xin — 0.065 (dotted) and fitted curve (solid) with AEc ~ 55 
meV. (Sample #298, T=300K). 
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FIG. 11: Position of bound energy level Ei in strained 
Ina;Gai_2, As/GaAs-QWs (triangles up for w < 7.2 nm and 
triangles down for w > 7.2 nm). Solid line — calculated de- 
pendence El of xin in the assumption w = 7.2 nm and for the 
doping concentration as in the sample #300. Dash-dot line 
— the corresponding position of Fermi level. 



remains significant to provide an excess of apparent car- 
rier concentration in QW region over the dopant value in 
all samples, except for x — 0.065. 

From Fig. it can be seen that there is no bound 
energy level for x < 4%. Indeed, the weak doping in the 
adjacent to QW spacers leads to an additional conduction 
band bending near the QW, which lifts the energy level 
up. The effect of disappearing bound level does not exist 
if the spacers are absent (at least, down to extremely low 
X, about 1%, when errors in numerical calculations begin 
to occur). 

In Fig. El the results on conduction band offsets in 
strained Inj,Gai_a;As/GaAs-QWs obtained during the 
numerical fitting to the experimental C-V characteris- 
tics are presented. The "recommended" curve from the 
above mentioned review^ is also depicted. The "rec- 
ommended" values of band offsets are higher by about 
25% in comparison to our results for strained quantum 
wells. The origin of this difference most probably is the 
presence of significant degree of elastic strain in pseudo- 
morphic InGaAs on GaAs. The deformation potential in 
compressively strained InGaAs modifies the energy line- 
ups of the heterostructure. Estimations on the base of 
model-solid theory'^'' predict an increase of the band gap 
of compressively strained InGaAs in comparison to the 
value AEg for bulk material (this effect again is of the 
order of 25%). So, the absolute values of band discon- 
tinuities are smaller in compressively strained quantum 
wells than in relaxed single heterostructures or thick dou- 
ble heterostructures. Recent results of other researchers 



EE c (eV) 
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FIG. 12: Conduction band offsets in strained 

Ina;Gai_j, As/GaAs QWs as a function of composition. 
Dashed line — as recommended in Ref. |3- Open circles — 
the latest results of capacitance and optical investigations on 
strained In^Gai_^As/GaAs-QWs: 1 - Ref . [H 2 - Ref. HI 
S-Ref.m 4-Ref.|il 



on strained QWs, mainly obtained by capacitance tech- 
niques (also depicted in Fig.lT^. are in reasonable agree- 
ment with ours, but exhibit significant scattering. 

The experimentally obtained dependence AEc — f{x) 
is close to a straight line with only little bowing. From 
fitting the curve to a parabola we propose the expression 
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AEc{x) = O.SUx - 0.21x2 fQj. ^j^g conduction band off- 
sets in strained In2;Gai_a;As/GaAs quantum wells in the 
composition range < x < 0.3. 

V. SUMMARY 

Aiming to get accurate and precise values for conduc- 
tion band offsets, a set of high quality samples containing 
strained Ina;Gai_a;As/GaAs quantum wells was grown in 
the composition range 0.06 x ^ 0.29. Specially for 
C-V-mcasurcments a constant impurity concentration in 
the cladding layers was maintained during the growth 
in order to eliminate uncertainties in subsequent numer- 



ical simulations. A fitting procedure of experimentally 
obtained apparent concentration profiles has been imple- 
mented using self-consistent solution of Schrodingcr and 
Poisson equations. All important information about the 
properties of the quantum well structures was derived: 
the majority carrier profiles, the positions of energy lev- 
els, corresponding wave functions, profile of the conduc- 
tion band bottom, as well as the dependencies of the 
above mentioned parameters on the applied electric field. 
The presence of only one bound level was discovered in 
all samples. The conduction band offsets in strained 
In2;Gai_2;As/GaAs quantum wells follow the expression 
AEcix) = 0.814X- 0.21x2. 
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